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The explanation is presented for the temperature dependence of the fourfold intensity peak split- 
ting found recently in diffuse scattering from the disordered CU3AU alloy. The wavevector and 
temperature dependence of the self-energy is identified as the origin of the observed behaviour. 
Two approaches for the calculation of the self-energy, the high-temperature expansion and the en- 
expansion, are proposed. Applied to the CU3AU alloy, both methods predict the increase of the 
splitting with temperature, in agreement with the experimental results. 
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Recently, in the first in situ experiment to resolve the 
fine structure of the equilibrium diffuse scattering inten- 
sity from the disordered CU3AU alloy, Reichert, Moss and 
Liang jjj have observed a marked temperature depen- 
dence of the fourfold splitting of the (110) short-range 
order (SRO) diffuse intensity peak. The separation of the 
split maxima changed reversibly, increasing with temper- 
ature. The same behaviour of the splitting was also found 
in H by analysing results of the Monte Carlo (MC) sim- 
ulations for the Cuo.856Alo.i44 alloy ||. The peak split- 
ting (Fig. [j]) is attributed to the indirect interaction of 
atoms via conduction electrons in an alloy whose Fermi 
surface has flat portions; the effective interatomic pair in- 
teraction itself has split minima in the reciprocal space, 
and their location is determined by the wavevector 2kp 
spanning these flat portions of the Fermi surface 0] . As 
indicated in [jjj , current theoretical approaches fail to ex- 
plain the observed behaviour. Indeed, the standard ap- 
proximation for the SRO diffuse intensity, the Krivoglaz- 
Clapp-Moss (KCM) formula j§, is 
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where /(k) is the intensity in Laue units, c the concen- 
tration, [3 = l/T, T the temperature in energy units and 
V(k) the Fourier transform of the combination Vij = 

(V4 A + V ij B )l 2 - V 4 B of potentials V t f with which an 
atom of type a at site i interacts with an atom of type 
(3 at site j. Eq. ([!]) predicts that positions of the J(k) 
peaks coincide with those of the corresponding minima 
of V(k); therefore, the splitting does not depend on T, if 
it is assumed that V^(k) is T-independent. This assump- 
tion is justified at least as far as positions of the V(k) 
minima are concerned, since the 2kp value should not 
change over the considered temperature range [Q. Be- 
sides, the MC calculations || in which the increase of the 
splitting with temperature was found |^] were carried out 
for the T-indepcndent pair interaction parameters. On 
the other hand, the cluster variation method || which in 
most cases leads to a significant improvement of the re- 
sults in comparison with the KCM approximation |t]], is 
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FIG. 1. Schematic reciprocal-space picture of scattering 
from the FCC alloys discussed in the text. Large circles repre- 
sent the Bragg reflections, while small ones (forming charac- 
teristic crosses) correspond to the split diffuse intensity peaks. 



practically inapplicable here, since interactions between 
atoms at distant lattice sites are involved (see below). 

The aim of the present Letter is to propose the theory 
of the temperature-dependent peak splitting observed in 
the CU3AU alloy. We begin by noting the exact expres- 
sion for the SRO diffuse scattering intensity 0], 
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where £(k) is the so-called self-energy which depends not 
only on k, but also on c and T. In the KCM approxi- 
mation, however, £ is k- and T-indepcndcnt, as follows 
from the comparison of Eqs. (Q) and (^): 
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Below we consider the J(k) profile along one of the lines 
containing split peaks, e.g., the (hlO) line, and concen- 
trate on two peaks around the (110) position. The peak 
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FIG. 2. Shift of the intensity peak position as a result of 
the wavevector dependence of the self-energy. The latter is as 
found for the CU3AU alloy (see text). 



positions kj (k is the deviation of the wavevector from 
the (110) position along the (hlO) line) are determined 
by the condition d k I = which gives 



2d k V = Td k Y, 
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Eq. ([!]) means that the k-dependence of £ leads to the 
shift 5k — kj — ky of the peak position with respect to 
the position kv of the corresponding minimum of V(k) 
(Fig. ||); kv is the solution of the equation d k V = 0. 
Furthermore, the right side of Eq. (0) is a function of T, 
while its left side is T-independent. The I(k) peaks will 
therefore change their positions with temperature. 

At sufficiently high temperatures the behaviour of the 
splitting can be analysed by using the high-temperature 
expansion (HTE, in powers of (3V) for S. The second- 
order HTE approximation gives || 
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Here £<j is the diagonal and £ 0( j the off-diagonal part 
of S. In this approximation Eq. (M) reduces to 



d k V=(l- 2c) 2 13 d k W , Wij = V2 
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The right side of Eq. (|6|) is small due to the prefactor 
/?, and its solution ki deviates little from kv- It is then 
sufficient to expand d k V and d k W in powers of the shift 
5k and retain only linear terms, 

d k V(k) = (d 2 k V) kv 5k , (d 2 V) kv > , (7a) 
d k W{k) = {3 k W) kv + {d 2 k W) kv 5k . (7b) 

Substituting Eqs. (f/1) into Eq. (||) and neglecting the last 
term in Eq. ( [7b| ) because of the smallness of the right 
side of Eq. (pi) lead to the result 



TABLE I. Coefficients B f and Cf (Eqs. (|l0|)) for the in- 
teraction V(k) and related function W(k) as calculated from 
the inverse MC interactions Vi m „ JlTJ for the CU3AU alloy at 
different temperatures. 
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Eq. (||) shows two scenarios for the temperature depen- 
dence of the splitting, depending on the sign of (d k W) kv . 
The first is the increase of the splitting with temperature 
discussed above. Apart from that, the theory predicts 
that the decrease of the splitting with increasing tem- 
perature is also possible; such temperature dependence 
has not yet been observed. The absolute value of 5k de- 
creases as T _1 with temperature. The shifts of the two 
peaks have opposite signs and the same absolute values, 
and the wavevector dependence of £ determines whether 
the splitting increases or decreases with temperature. In 
the case of an equiatomic alloy (c = 0.5) the second- 
order contribution to d k T, is zero, and the temperature 
behaviour of the shift is defined by the higher-order terms 
in the HTE for the self-energy. 

We now apply the HTE to the CU3AU alloy for which 
sets of first 8 inverse MC interactions Vi mn were ob- 
tained at different temperatures using SRO parameters 
o-imn available in the literature [ |l(| . Despite its extended 
range, the behaviour of the interaction along the (hlO) 
line is simple; the Fourier transform /(k) of an arbitrary 
FCC matrix /y with non-zero elements for the first 20 
coordination shells Imn has along this line the form 

f(k) = Af + 2Bf cos27rfc + 2(7/ cos47rfc , (9) 

where the relevant coefficients Bf and Cf are 

Bf = /200 — 4/211 + 4/220 + 4/222 — 8/321 

+4/420 - 4/ 332 + 8/422 - 8/521 + 4/442 , (10a) 

Cf = /400 — 4/411 + 4/420 + 4/422 — 8/431 

+4/440-4/433 + 8/442, (10b) 

and k is measured in the reciprocal lattice units (r.l.u.). 
These coefficients for functions V{k) and W(k) are shown 
in Table |. All four quantities scatter widely, which is the 
result of relatively low accuracy of the inverse MC inter- 
actions for the CU3AU alloy discussed in |l(|. In par- 
ticular, the accuracy and/or number of interactions are 
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TABLE II. AE coefficients B f and C f (Eqs. @) for the 
self-energy E(fc) as calculated within the 10-shell approxima- 
tion (except for set 16; see text) from the experimental sets 
of the SRO parameters for the C113AU alloy at different tem- 
peratures. 
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Ref. 
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2 
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10 
11 
12 
13 
14 
15 
16 
17 



669 
693 
748 
958 

1023 
669 
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0.1742 
0.0905 
0.0915 
0.0268 
0.0231 
0.2172 
0.0787 
0.0592 
0.0130 
0.0120 
0.2086 
0.1372 
0.0530 
0.4137 
0.2053 
0.1791 
0.0911 



0.0079 
0.0119 
0.0013 
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0.0006 
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0.0004 
0.0006 
0.0144 
0.0037 
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insufficient for the correct description of the split mini- 
mum of V(k). The splitting in V(k) occurs when Cy > 
and I By I < 4Cy, and this is so only for sets 1 and 3. In 
cases 4 and 6 Cy — 0, since the eighth interaction V400 
(the only one contributing to Cy) was not included in 
the corresponding sets. In addition, for cases 6 and 8 
By > 0, so that the split minimum would occur around 
the (|10) rather than (110) position. Nevertheless, de- 
spite low accuracy it is seen that Bw is positive and Cw 
is non-negative. The only exception is set 6, where Bw 
acquires very small negative value. However, the second 
set of interactions (set 7, with non-zero V400) obtained us- 
ing the same SRO parameters leads to positive values of 
both Bjy and Cw- In the case Byy > 0, Cw > function 
W(k) has a maximum at the (110) position which is much 
wider than the magnitude of the peak splitting; \ki\ val- 
ues observed in |j| were quite small (less than 0.1 r.l.u.). 
As a result, the derivative dkW is positive for the left and 
negative for the right minimum of V(k) (Fig. ||), and at 
any finite temperature the intensity peaks are shifted to- 
wards the (110) position. The absolute value of the shift 
increases with decreasing temperature, so that the split- 
ting increases with temperature. 

The applicability of the HTE, similarly to that of the 
KCM approximation, is limited to the case of sufficiently 
high temperatures. To deal with moderate temperatures, 
we introduce here another approach which leads to the k- 
and T-dependence of S, its expansion in powers of SRO 
parameters ctij (hereafter the a-expansion, or AE). Two 
non-zero orders of the AE for £ <2 were calculated |8| in 
the framework of the 7-expansion method (GEM) 
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The expression for then comes from the sum rule 
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(here the integration is carried out over the Brillouin zone 
of volume to), Eq. @ being one of the AE (or GEM) 
equations: 
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Note that the sum of the first two terms corresponds to 
the spherical model (SM) for SRO |l2| , which is the zero- 
order approximation for the AE and GEM; in the SM the 
self-energy is diagonal (a = b = 0). The difference be- 
tween the AE and GEM lies in the choice of the expansion 
parameter (s). The GEM parameter is 7 = exp(— 1/£), £ 
being the dimcnsionlcss correlation length, and terms in 
the diagrammatic expansion for the self-energy are se- 
lected according to the total length of all lines in the dia- 
grams, where the line connecting sites i and j represents 
ctij ■ The GEM is based on the assumption that the cor- 
relations decrease rapidly with distance; this assumption 
is invalid here because distant interactions are essential. 
The AE uses ctij themselves as the expansion parameters; 
the terms are chosen according to the number of lines in 
the diagrams (i.e., the powers of ctij), since all ctij are 
sufficiently small. The GEM was successfully applied to 
both the direct and inverse problems of alloy diffuse scat- 
tering [§[l3) , leading to reliable results everywhere except 
in the vicinity of the instability point. Based on our ex- 
perience with GEM, we can as well expect the AE to be 
quite accurate at almost all temperatures. 

Applying the AE to the CU3AU alloy, we calculate 
coefficients Bs and Cs combining Eqs. (10) and jll| ) 
and using available sets of experimental SRO parame- 
ters @-(L§. Their values for the case of first 10 shells 
included in the AE approximation ( |j"T| ) for the self-energy 
are given in Table [n| (5-shell AE approximation was used 
for set 16, since only 5 SRO parameters were determined 
in p8|). Inclusion of additional shells does not alter the 
results. In all cases Bs is positive and Cy. is non-negative, 
so that, as before, S(fc) has a maximum at the (110) posi- 
tion which is very wide in comparison with the peak split- 
ting, and the intensity peaks are shifted towards this po- 
sition. Contrary to the case of the HTE, the explicit tem- 
perature dependence of the AE self-energy is unknown, 
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FIG. 3. Values of g (Eq. (H)) vs. temperature calculated 
using the data from Table |j Jl4|] - dots, solid line; [jl5| - 
squares, dashed line; Jlq] - upward triangles, dot-dashed line; 
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downward triangles, dotted line; |18| 
rightward triangle. 



leftward triangle; 



since the SRO parameters in Eqs. (|11| ) are complicated 
functions of temperature. To find the temperature be- 
haviour of the splitting, we use the data from Table |l| 
and plot in Fig. || against temperature a quantity 
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which characterises the temperature dependence of the 
right side of Eq. ([!]) at small k. It is seen that g is a de- 
creasing function of T , which corresponds to the increase 
of the splitting with temperature. Its temperature de- 
pendence is particularly strong in the range below 800 K, 
where the intensity profile was measured in pj]. 

To summarize, we have presented the explanation of 
the temperature dependence of the Fermi surface-induced 
diffuse intensity peak splitting found recently for the 
CU3AU alloy. The wavevector and temperature depen- 
dence of the self-energy is understood to be the origin of 
this behaviour. The proposed theory is able to describe 
the observed increase of the peak separation with tem- 
perature; in addition, it also predicts the possibility for 
the splitting to decrease as temperature increases, the 
behaviour which has not yet been found. Two methods 
for the calculation of the k- and T-dependent self-energy, 
the HTE and the AE, have been proposed. Applied to 
the existing experimental and inverse MC data for the 
CU3AU alloy, both methods predict the increase of the 
splitting with temperature, in agreement with the exper- 
imental findings. However, the HTE is not expected to 
be reliable when applied to alloys at realistic tempera- 
tures, so that the AE approach is preferable. Despite the 
seeming complexity of the problem (the interaction in- 
volves many coordination shells), the theoretical analysis 
proves to be surprisingly simple. 
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